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Abstract 

We compute the exact S-matrix and give the Bethe ansatz solution 
for three sigma-models which arise as subsectors of string theory in 
AdSs x 5*^: Landau-Lifshitz model (non-relativistic sigma- model on S*^), 
Alday-Arutyunov-Frolov model (fermionic sigma-model with su(l|l) 
symmetry), and Faddeev-Reshetikhin model (string sigma-model on 
x R). 
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1 Introduction 

According to the AdS/CFT correspondence [1^ solving four-dimensional Åf = A super- 
Yang-Mills (SYM) theory amounts to quantization of type IIB superstrings on the 
AdSs x 5*^ background. This potentially easier problem has so far resisted solution, 
partly because the sigma-model on AdSs x S*^ is of the Green- Schwarz type and has 
well-known difficulties in the conformal gauge. The sigma model, however, is integrable 
[3] which gives us hope to solve it with the help of Bethe ansatz. There is a mounting 
evidence that the non-perturbative spectrum of AdS/CFT is indeed described by some 
sort of Bethe equations [HEI- The strengest evidence comes from studying the spec- 
trum of anomalous dimensions in the SYM theory jEl 13 E] and from analyzing classical 
Solutions in the sigma-model [31101111], which describe spinning strings in AdSs x 

The Bethe ansatz E] is a common method to solve integrable models, which 

encodes the whole spectrum of the system in a set of algebraic equations. The Bethe 
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equations can be interpreted in terms of factorized particle scattering, where the "Bethe 
particles" do not necessarily coincide with the physical degrees of freedom of the theory. 
In such an interpretation, the Bethe equations arise as quantization conditions for the 
momenta of the particles in a box of hnear size L: 

where A{pj,pk) is the phase shift experienced by the jth particle when it scatters off the 
kth particle. The energy of a Bethe state is the sum of the single-particle energies: 

j 

The energy spectrum is thus completely determined by the two-body scattering and the 
one-body dispersion relation. There are no genuine multi-body interactions [TH] . 

The SYM theory is an example where the Bethe particles are very different from 
the excitations in the four-dimensional space-time. Here, the Bethe equations, which 
determine the planar spectrum of anomalous dimension, were derived by interpreting 
single trace operators in SYM as quantum states of abstract spin chains [Hj. The Bethe 
particles are waves that propagate along the spin-chain associated with a local operator 
in the SYM, but more recent results suggest that the story might be different (and 
far more interesting!) at the non-perturbative level. On the string-theory side, the 
evidence in favor of the Bethe-ansatz structure of the spectrum comes from the two 
sources: (i) The classical solutions of the string sigma-model can be parameterized by 
the integral equations of Bethe type [Hl CD] and (ii) the leading quantum corrections 
in the near- BMN limit of the AdSs x S'^ geometry (THj can be parameterized by a set 
Bethe equations [121 1201 • The quantum string Bethe equations, which were conjectured 
on the basis of these two observations 1201 IHl 1^ , receive corrections at higher orders 
in the sigma-model coupling [22\ and thus contain infinitely many unknown parameters. 
Deriving Bethe equations for quantum string in AdSs x from first principles is an 
open problem. In particular it is not quite clear what degrees of freedom of the string 
are represented by the Bethe particles. 

The string sigma-model on AdSs x 5*^ is a rather complicated two-dimensional field 
theory. Similar but simpler systems, such as the Osp(2m + 2|2m) coset sigma-model, 
were solved by Bethe ansatz j2Sl, because many two-dimensional factorized S-matrices 
are known exactly [Mj. For other models a direct relationship to spin chains can be 
established 24 . The idea that we would like to put forward in this paper, and test 
on a number of simplified models, is to derive the Bethe equations for quantum strings 
by explicitly comp uting the two-body scattering matrix on the world-sheet. We use 
relatively simple and rather standard methods of quantum field theory to do that. This 
will allow us to obtain the Bethe equations for several two-dimensional field theories 
that arise as reductions of the AdSs x 5*^ sigma-model. We should mention that the 
S-matrix approach was quite successfully used on the gauge-theory side of the AdS / CFT 
correspondence [201 123 1201 ■ 
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We demand (rather than prove) quantum integrability and factorization of the S- 
matrix as necessary prerequisites. With these assumptions in mind, our method provides 
the means to derive the Bethe equations of an integrable theory at the quantum level. A 
more rigorous approach to Bethe ansatz is the quantum inverse scattering method j2Z] 
which fully exploits the rich algebraic structure associated with integrabihty. 

The first model we consider is a non-relativistic sigma model on 5*^ (the Landau- 
Lifshitz model). This model describes fast moving strings on x R and arises as a 
low-energy effective theory of the Heisenberg ferromagnet (the connection that plays an 
important role in the AdS/CFT correspondence |2H1)- We derive the Bethe equations 
for the LL model in Sec. 01 by direct computation of the S-matrix. In Sec. HJwe consider 
the fermionic sigma-model that arises from the su(l|l) reduction of the AdS string [29]. 
Finally, in Sec. Øwe consider again string theory on x R, but this time without making 
the low energy approximation. This model was introduced by Faddeev and Reshetikhin 
and is closely related to the su(2) principal chiral field, whose Bethe-ansatz solution 
was obtained in jni]. The lattice-regularized quantum version of the FR model was 
solved in I^^OJ. We derive the Bethe equations directly in the continuum. 



2 Preliminaries 

The general idea is to derive the Bethe equations for quantum strings from scattering 
computations in the world-sheet theory. Let us sketch this idea for the example of a 
single scalar field on the world-sheet. In the usual Hamiltonian approach to the Bethe 
equations, one constructs eigenstates which have the structure of scattering plane-waves. 
In the two-particle sector such an eigenstate of the Hamiltonian looks like 

\pp') = f dxdx' [e{x' -x) + 9{x - x')S{p,p')] e^P"+*P'"' ip\x)ip\x') |0) . (2.1) 

J ^ v ' 

It is parameterized by two momenta p and p'. We will always label the momenta such 
that p > p' , so that the first term in the wave function is the incoming wave, and 
the second term is the scattered wave. Imposing the periodicity condition on the wave 
function: x(0,a;') = x{L,x') and x(^;0) = x(z, L), one finds that the momenta of the 
particles are quantized according to 

e'^P = S{p',p) , e^^^P' = S{p,p'), (2.2) 

where S{p',p) = 1/S{p,p'). This is a particular case of (jl.lll . for which integrability is 
actually not required. The distinguishing feature of integrable models is that the multi- 
body wave function is two-particle reducible jnSIISl- The periodicity condition for an 
arbitrary multi-particle eigenstate of the Hamiltonian is then expressed in terms of the 
two-body phase shifts as in 

Diagonalization of the Hamiltonian, however, is not the most efficient way to calculate 
the scattering phase shifts in field theory. It is much easier to compute S{p,p') as 
the matrix element of the infinite time evolution operator S between the two-particle 
scattering states: 

{kk'\S\pp') = S{p,p')5MXk') . (2.3) 
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Here we have introduced the notation 

5^{p,p\ k, k') = (27r)2(5(p - k)S{p' - k') ± 5{p - k')5{p' - k)) . (2.4) 

This factor represents the conservation of individual momenta during the scattering 
process. In two dimensions, energy and momentum conservation allows two particles 
only to exchange their momenta. The relative sign between the two terms is plus for 
bosons and minus for fermions. 

We will compute the S-matrix using Feynman diagrams. Note that the usual Feynman 
rules calculate the matrix element A^(p, p', fc, fc') as defined in 

{kk'\{S - l)\pp') = iM{p,p ,k,k') {2Tif 5^^\p^ + p"" - k^ - k"") . (2.5) 

The energy-momentum conserving delta-function is different from ()2.4|) by a Jacobian 
l/{de/dp — de/dp'), which has to be tåken into account when extracting the phase shift 
from the diagrammatic calculations. 

In what follows we make use of ()2.3|1 and derive the Bethe equations for three different 
theories that are of potential relevance for quantum strings in AdSs x S^. 



3 Landau-Lifshitz model 

The LL model is defined by the action 



S 



d x 



Ct{n) - -{d^nf 



where n is a three-dimensional unit vector: 



n 



1 . 



The action contains the non-local Wess-Zumino term 



di Sijk rii d^rij dgUk 



(3.1) 



(3.2) 



(3.3) 







and is of the first order in time derivatives. The equations of motion that follow from 
dSH) are 

dtUi = Eijk Tij dlrik . (3.4) 

The LL equation is completely integrable ^2^. The quantum inverse scattering method 
for the LL model was discussed in |34] . 

The LL model is the low-energy effective field theory of the Heisenberg ferromagnet 
with Hamiltonian 

^ a.+i) , (3.5) 



r 



1=1 



which in the AdS / CFT context arises as the one-loop mixing matrix of scalar composite 
operators ii{Z^~^W^^ + permutations) The Heisenberg equations of motion for 
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reduce to ()3.4j) in the continuum limit L — oo if the spin operators are formally 
replaced by unit c-number vectors: ?;(r) n{T,a) with a = 2ttI/L. Alternatively, the 
action in the coherent-state path integral of the Heisenberg model in the continuum limit 
becomes 



S 



L 
2^ 



dr 



2-E 



do 



C{n) 



A 

8L2 



(3.6) 
]nH/X. The 



This is the same as 1)3.11) after the following rescalings: a 
spacial coordinate x in ()3.1)1 now has periodicity L. 

The LL model can be also derived from classical string theory on S^xR in the limit of 
fast-moving strings. The details of the derivation, together with the precise matching to 
the mixing matrix for long operators in the SYM, can be found in |2S|. Here we would like 
to view the LL model as a (l+l)-dimensional quantum field theory. We derive the Bethe 
equations for the non-perturbative spectrum of the LL model with the help of relatively 
simple perturbative calculations, which are similar to the perturbative calculation of 
the S-matrix in the closely non-linear Schrodinger model [33]. Quantum-mechanical 
(Hamiltonian) perturbation theory for the LL model was developed in [321 EZI- Here, we 
will use Feynman diagrams which greatly facilitates the calculation of the S-matrix. 

In order to develop perturbation theory we first re-write the WZ term in the local 
form. For that we will need some properties of the WZ action, which we review here 
following [3H]. The WZ term in ()3.1|) can be written as 

WZ\n\ := f dtCtin) = 



Eijk riidnj A duk 



where dn = dtndt + dgndO. A short calculation shows that 



-'ijk 



riidrij A dnf, 



drii A dn2 



^ / nidn2 — n2dni 
V 1 + ^3 



(3.7) 



(3.8) 



2 ^ ^ ns 

This identity allows one to write the WZ action in a local form, which comes at the price 
of losing manifest SO (3) invariance: 



WZln] 



- / dt 



nin2 - ^2^1 



(3.9) 

1 + ^3 

The next step is to sol ve the constraint ()3.2|) by expressing in term of rii and n2- 
It is also convenient to make a field redefinition which gets rid of the non-linearities in 
the kinetic term [SH], and also to combine ni, n2 into a single complex scalar, since the 
complex field has the canonical non-relativistic propagator. The significance of this fact 
will become clear later. So we define 

ni + in2 



2m 



V2 + 2^3 

Performing this change of variables in ()3.9j) . 1)3.1)) . we find: 



(3.10) 



S 



d x 



1 - 



(3.11) 
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Figure 1: Pole prescription in Landau-Lifshitz model. 



This action describes an interacting field theory of a single scalar field. Its non-relativistic 
character leads to some important non-renormalization properties. 
First of all, the ground state is annihilated by the field operator: 

ip{t, x)\0) = . (3.12) 

Since the equations of motion are of the first order in time derivatives, the field operator 
in the interaction picture is expanded in negative-frequency modes only: 

^it,x) = l^a,e-^^''^^^- ' V'*(^>^) = /^«^e^"'*~^'^ (3-13) 
where ap, a]p create and annihilate a particle with momentum p and energy^ 

e{p) = . (3.14) 

The operators ap and aj, obey canonical commutation relation normalized as 

[ap,al]= 271 6{p-p') . (3.15) 

Since the ground state is annihilated by ip{t,x), the particles do not travel back- 
wards in time and the propagator, accordingly, has only one pole in the momentum 
representation: 

L>(t,x) = (0|T(^(t,a;)^*(0,0)|0) = . < . 



(3.16) 

(27r)^ uj — p^ + ie 

The pole prescription, see Fig. or again the fact that ^{t,x) annihilates the vacuum 
makes the coordinate-space propagator purely retarded: 

D{t,x)=e{t)J-e- . (3.17) 
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Figure 2: Non-renormalization theorem. Any closed loop of likewise oriented propagators 

Dmn = Ditm - tn,Xm - Xn) VailisheS. 



Å 









x x 


XXX 



Figure 3: Poles for a closed loop. All propagators in a closed have their poles in the lower 
half plane. Hence the integral (|3.19|) over the energy flowing around the loop vanishes. 



We can now prove the following "non-renormalization theorem": any diagram that 
contains a closed loop with arrows in the same direction vanishes, cf. Fig. |21 In the 
coordinate representation this follows from the fact that at least one propagator in the 
loop has a negative time argument. Consequently, 

D{t2 - ti, X2 - Xi)D{t^ - t2, X3 - X2) ■ ■ ■D{ti - t„, Xi - Xn) = . (3.18) 

In the momentum space representation, the integrand has poles only in the lower half- 
plane of complex u. The integral is then zero by the contour argument, cf. Fig. El 



/duj i i i f3 19) 

This theorem has three important consequences: 

• The ground state energy is not renormalized: i?vac = 0. This is consistent with 
the fact that the ferromagnetic vacuum is the exact zero-energy eigenstate of the 
Heisenberg Hamiltonian ()3.5|) (cf. the discussion in 



^Recall that we have rescaled the time variable by Stt^/A. Hence, the energy is the frequency 
multiplied by A/Stt^ 
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Figure 4: Generic loop diagram for the two-body S-matrix. 




Figure 5: Cutting a generic diagram. Due to charge conservation the number of future 
directed propagators minus the number of past directed propagators has to be the same at any 
moment in time. Since any past directed propagator is identicaUy equal to zero, the number of 
propagators at any cut is the same and is equal to the number of external incoming/outgoing 
legs. 

• The one particle Green's function is not renormahzed. Hence, the dispersion rela- 
tion ()H.14j) does not receive quantum corrections. 

• The two-body S-matrix is given by the sum of bubble diagrams (Fig. HJ. 

These properties are almost obvious. A formal proof can be given by cutting a generic 
diagram and counting intermediate propagators, as illustrated in Fig. |S1 

The fact that the two-body S-matrix is determined by the sum of bubble diagrams in 
Fig. 0] has far reaching consequences. Since these diagrams contain only quartic vertices 
we may truncate the non-polynomial action (jH.llll at the fourth order in the fields^: 

L='- {^*dt^ - dt^*^) - \d^^f - I {{^^d^^f + (4<^V)'] + . (3.20) 

This is a very important simplification, which makes the all-loop computation feasible 
and allows to forget about other non-linear terms in the action. In ()3.20j) we have 
introduced a formal expansion parameter g to make the power-counting of perturbative 
series more transparent. We will set 5» = 1 at the end of the calculation. In fact, 
observable quantities cannot depend on (7, since this parameter can be eliminated by 
rescaling t, x and (/?. 

Let us now derive the Bethe equations for the LL model ()3.ip respectively ()3.11|) . i.e. 

^The resulting Lagrangian is very similar to the one of the non-hnear Schrodinger model given by 
£ = ^ {ip* dtif — dt^* Lp) — \dxpf' — g{<f*(p)'^. For comparison we recall the S-matrix of this model: 
S{p,p') — pZpi^lg 1 which can be computed by the same technique as we use here |85| . 
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compute the two-particle S-matrix. The two relevant vertices, written in ()3.2()j) . are 

A; v 




^ ^ =2zghk[. 



(3.21) 
(3.22) 



(3.23) 



(3.24) 



The two-particle in- and out-states are defined as 

\pp') = a\p)a\p')\0) , {kk'\ = {0\a{k')a{k) . 
The non-scattering part is given by 

{kk'\pp')=54p,p',k,k'), 
where 5+ was defined in fl2.4|) . At tree-level we need to evaluate the following expression 

{kk'\S\pp')\^ = {kk'\ (^-|^ Jdtdx[{ip*d,ipf + {d,ip*ipy]\pp') . (3.25) 

There are four identical Wick contractions between the states and each vertex leading to 
the factor 4{—kk' — pp'). The integration over t and x imposes energy and momentum 
conservation, which can be written as 

{2TTy6ip^+p'''-e-k''')5ip + p' -k-k') = —L—5+{p,p\k,k') . (3.26) 

2{p-p') 

Multiplying all factors together we find 

Tin' 

(3.27) 



{kk'\S\pp') =2zg^6+{p,p',k,k'). 

9 P-P 



Including the non-scattering part ()3.24|1 . we can read off the S-matrix from comparison 
with ()2.3|1 and find to first order 



S{p,p') = l + 2zg^ + 0{g 



(3.28) 



p — p 

In App. E]we compute all higher order corrections and find for arbitrary n, cf. flA.l 



{kk'\S\pp') =2{igy 

9" 



pp 



p — p 



S+{p,p',k,k') . 



(3.29) 



Summing this result for all n and adding the non-scattering part yields the exact S-matrix 
for the LL model 



s{p,p') = i+2j2{^9r 



n=l 



pp' 



I \ n 



p — pl 



i _ _L 

p p' 



(3.30) 



where the expansion parameter g can now be set equal to unity. 

The Bethe equations resulting from this S-matrix are most conveniently written in 
terms of the spectral variable u = l/p: 



Uj — Uk — I 



n ^' ■ (3.31) 

X X 7/ ■ 111 7 

The energy and the momentum are then given by 



E 



J i 



The Bethe equations for the spectrum of the Heisenberg Hamiltonian ()3.5|) are very 
similar ^1]. They are obtained by replacing l/u ^ ir— 2 arctanu and 1/u'^ —>■ 1/(m^+1/4) 
in the Bethe equations and the dispersion relation above. The difference disappears at 
large u (small momenta), as expected. In particular, the spectrum of the low-energy 
modes with p~l/n~l/L, is the same for both models and can be calculated from the 
approximate equation 

L 2 

27m,- = V , (3.33) 

«i t^j - Uk 

up to and including finite-size 0{1/L) corrections [Q. As shown in jSTj, the spectra for 
the two models start to deviate from one another at order 0{1/L'^). 

The classical limit is achieved by assigning a macroscopic number of Bethe roots to 
a finite set of mode numbers n/. Then fl3.33|l becomes an integral equation jSH]: 

/^ = l-2vrn„ (3.34) 
J x — y x 

where 

i 

It can be shown that ()3.34|) describes all time-periodic classical solutions of the LL 
equation ()3.4j) [9 . This classical limit of the Bethe equations is obviously the same for 
the LL and the Heisenberg models. 

The spectrum that foUows from Bethe equations (j3.3Hl indicates a rather unexpected 
instability. The bound states of elementary excitations are described in integrable field 
theories by string solutions of the Bethe equations. For instance, the 2-string configura- 
tion is 

i 

«1,2 = v±-. (3.36) 

It becomes an exact solution of the Bethe equations in the strict thermodynamic {L 
oo) limit. The energy and the momentum of the 2-string are: 



A 



1 1 



1 2v 



A v'-\ 



2 ' 



2-string , T ■ 9,1 

V + ^ V - ^ + ^ 



(3.37) 
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Figure 6: Analytic structure of the two-particle Green's function ()3.40|) . 



The momentum of the 2-string is always smaller that 2, so the 2-string solutions of the 
Bethe equations describe two kinds of excitations with the dispersion relations 

e±{p) = ^ (±2v/4-p2_4 + p2j . (3.38) 

The two branches arise from v ^ 1/2. The £_ branch has negative energy. 

The existence of negative-energy states seems to signal a vacuum instability. Such 
an instability is absent in the Heisenberg spin chain, and its appearance in the LL model 
is rather puzzling. We have no simple explanation for this phenomenon, but we can 
demonstrate that the negative-energy states do arise as poles of the Green's functions. 
Consider, for instance, the Green's function of a composite operator 

O = {v>*d,if'f + {d,^*^f , (3.39) 

which is computed in App. IA. 31 



(0| T 0{t, x)O{0, 0) I 0) = - / e-*^*+^^^ . (3.40) 

As a function of Q, it has a two-particle cut from the threshold at fl = P^/2 to infinity 
and two poles at g^fi = £±{P) (Fig. 



4 Alday-Arutyunov-Frolov model 

In this section we derive the Bethe equations for an integrable system of two-dimensional 
fermions which was recently found by Alday, Arutyunov and Frolov [22] ■ This model 
is a new integrable quantum field theory whose Bethe-ansatz solution is not known. 
Classical integrability of the model was demonstrated in [SH] by construction of a Lax 
representation for the classical equations of motion. The model arises as a consistent 
truncation of the classical AdSs x superstring to the su(l|l) subsector in the temporal 
gauge^. 

■^In the light-conc gauge the su(l|l) sector is described by free world-sheet fermions |40| . 
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The AAF model is a theory of an interacting Dirac (two-component complex) fermion 
in two dimensions"^ : 



27rL/^A 

S = — — / da I dr 

zn 



2 



° -y'^i^l;^a^|J^lJ^'^p^|J-^^^P^|J^pi>^^^P) ^^"^^ 

The action is manifestly Lorentz-invariant. The Lorentz metric is tåken in the (+1, — 1) 
signature and the notation x refers to both components x^) = (r, a). Below, we will 
use the scalar product a - b = g^^aabp and the "vector" product a x b = e'^^aabi3, where 
the epsilon tensor satisfies e^^ = eio = +1. The exphcit form of the Dirac matrices is 
given in ()B.2jl and the conjugate spinor is defined as ip = ip^j^. 

In the action ()4.ip all quantities {ip, a, r, L, A) are dimensionless. However, we prefer 
to write it in a way which allows us to assign canonical mass dimensions to the field and 
the coupling constants. Therefore we rescale the world-sheet coordinates by — > -j^x" 
and introduce the parameters 

2vr TT , , 

m = ^ , g = — ^ . (4.2) 

Then the space integral runs from to L and the Lagrangian can be written as^ 

C = — — m)ip 

For the purpose of power counting, we may assign the foUowing mass dimensions: [ip] = |, 
[L] = [x] = —1, [m] = 1, [g] = 0. Clearly, m plays the role of a mass parameter 
and g the role of a coupling constant. The model is integrable for any m and g. It 
is not renormalizable by power-counting, but as conjectured in [201, the symmetries 
(integrability being one of them) can take care of renormalization and render the model 
renormalizable in perturbation theory, the same way symmetries make two-dimensional 
non-linear sigma-models renormalizable. In this sense, the AAF model is a fermionic 
counterpart of the non-linear sig'ma-models. 

^We use the Lagrangian (5.7) from [5^, take the scaHng (5.6) into account and denote the angular 
momentum of the string by L instead of J. In addition we shifted a ^ a + l-KLj^fX to get positive 
hmits for the integration. Finally we drop the constant term in the action, which amounts to shifting 
the origin of the energy scale. 

^Now this theory resembles the massive Thirring model: £ = -0(1^ — m)i/; — ^ 'tp^°'tp 'ip'ya'4'. For 

taj^jj — 8' 

later comparison we recall the S-matrix of this theory: S{9,9') — — jrp-- We checked that this 

S-matrix can be reproduced by the same diagrammatic calculation as we do below for the AAF model. 



12 



One cannot compute the physical S-matrix of the AAF model directly by resumming 
diagrams, as we did for the LL model, because the usual relativistic propagator is not 
purely retarded and loop corrections do not cancel. This difficulty is not fatal and can 
be overcome by the trick that dates back to [H] and was used by Bergknoff and Thacker 
to solve the massive Thirring model P2]- The idea is to use pseudo- vacuum (instead 
of the true ground state) as a reference state for quantization. The pseudo-vacuum, by 
definition, is a state annihilated by the field operator: 

i;{x)\0) = 0. (4.4) 

All anti-particle levels in the pseudo-vacuum are left empty. Quantizing in the pseudo- 
vacuum can thus be interpreted as applying an infinite negative chemical potential to the 
system. The scattering S-matrix of excitations over the pseudo-vacuum can be computed 
exactly by essentially the same method as in any non-relativistic theory. This "bare" 
S-matrix can be then used to write down the Bethe equations. By solving the Bethe 
equations one can fill back the Dirac sea and reconstruct the true ground state. This is 
a pretty standard way to solve relativistic integrable field theories ^^l- In principle it 
allows one to derive the spectrum of physical states and their S-matrix [l^ from first 
principles. The price to pay is that the solution of the Bethe equations that describes the 
ground state and the solutions that describe physical excitations are already rather com- 
plicated. It is obvious that the filling of the Dirac sea drastically changes the spectrum 
of excitations and their S-matrix. 

By solving the free equations of motion we find the following mode expansion of the 
fields 

^{x) = I ^ [a\p,)uip,)e'P^-' + b\-p^)v{-p,)e-'P^-'] , (4.6) 

where po = ^/pl + m?- The spinors u{pi) and v{pi) are defined in App. IB. 11 The 
oscillators obey the standard commutation relation: 

{a{p,),a\p\)} = 2n5{p,-p\) , {b{p{),b\p',)} = 2t: 5{p, - p\) . (4.7) 

The operators a^{pi) and h'^{pi) create excitations with momentum pi and with energies 
-\-^/p{^^rn? and —\/pi + rn^- We define the pseudo-vacuum |0) of the theory as the 
state satisfying 

a(pi)|0) = 6(j9i)|0) = for all pi . (4.8) 

We call the excitations created by and h"^ on top of this state pseudo-particles with 
positive and negative energy, respectively. The physical vacuum \VL) is obtained from |0) 
by exciting all negative energy modes, i.e. by filling all holes in the Dirac sea. 

The major advantage of quantizing the theory in the pseudo-vacuum is the absence of 
anti-particles, which implies the same non-renormalization theorems as in the LL model 
of Sec. El Since the field operator %Ij[x) annihilates the pseudo-vacuum, the propagator 
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Figure 7: Pole prescription in Alday-Arutyunov-Frolov model. 



is purely retarded: 



Dix 



x 



(o|rv^(f)V^(f )|o) 



+ m) 



-ip-ix—x') 



(4.9) 



(27r)2 f 



■w? + ie(po) 



where e(po) = sign(po)e- The fact that D{x — x') vanishes for x° < can easily be 
seen from the pole prescription which is depicted in Fig. [7| We mention that this pole 
prescription is relativistic invariant because sign(po) is invariant under orthochronous 
Lorentz transformations. Now, as both poles are in the lower half plane, we conclude by 
precisely the same reasoning as in the previous section that the closed loop of likewise 
oriented propagators vanishes, cf. Fig. |21 The most important consequence of this result 
is the fact that the two-body S-matrix is entirely determined by the four-valent vertices 
alone. These are given by 



x p 1 (g) 7^ 



X= +^k' x Å^l 

A. p 



(4.10) 
(4.11) 



The unit matrix 1 (represented by the black dot) connects the p and k legs, whereas 
7^ (represented by the white dot) connects the p' and k' legs. The six-valent vertex is 
crucial for the factorization of the S-matrix and hence for the integrability of the theory, 
but it is of no relevance for the scattering of two pseudo-particles! 

In the foUowing we compute the S-matrix between two particles of type a"'' (which 
have positive energy). The scattering of particle of type 6^ can be obtained by analytic 
continuation to complex rapidities 6, defined by 



Po = m cosh 6 , Pl = m sinh 6 . 



(4.12) 



Real rapidities 6 = a & R parameterize the mass-shell of the positive energy particles 



and complex rapidities 6 



ITT 



a with a E R parameterize the mass-shell of the 
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Figure 8: Rapidity plane. It is convenient to work with a complex rapidity 9 as particles 
(lm ^ = 0) and anti-particles (lm 9 = ir) can be treated at the same time. 



negative energy particles 6^, cf. Fig. |H1 As a meromorphic function of the rapidities, the 
S-matrix describes the scattering of both types of particle at the same time. Hence it is 
sufficient to consider the following in- and out-states: 

\piP'i) = a\pi)a\p[)\0) , {hk[\ = {0\a{k[)a{h) . (4.13) 

The non-scattering part is given by 

{k^k[\p,p[)=S4p^,p[,h,k[), (4.14) 

where (5_ was defined in fl2.4|) . At tree-level we need to evaluate the following expression 



{k, k[\S\pip[)\^ = {kl k[\ l^-^j j c^2a;£"^[^a«#7'5/3^ - do^i^^pOfs^Pj^^P] \pip[) . 

(4.15) 

As the vertices do not have any symmetry, connecting the external lines with the vertex 
leads to 2^ = 8 different terms. One of them is given by 



IV? 



{p x p) {u{k[) ® u{ki)) (1 ® 7^) {u{pi) ® u{p[)) . (4.16) 



The others are obtained from this by exchange of the in-going particles {p p'), by 
exchange of the out-going particles {k k') and by choosing the second vertex {p' xp ^ 
k' x k). All changes are accompanied by a change of the overall sign. The spacetime 
integral imposes energy and momentum conservation, which can be written as 

(27r)'5{p + p'-k-k') = -É^_^S4pi,p[,k„k[) . (4.17) 

p' x p 

The 5-functions in S^{pi,p[, ki, k[), cf. ()2.4|) . can be used to convert the ^'s into p's. After 
this, the combination with the opposite relative sign S-{pi,p[, ki, k[) forms. Eventually, 
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the whole expression boils down to 

2iq 

{kl k[\S\pip[) =^{p' xp) (u{p[) u{pi)) (1 ® 7^ - 7=^ ® 1) (u{pi) ® u{p[)) 



p' X p 

(4.18) 



The scalar product of the spinors is given by (for useful identities see App. IB.1|) 

3 

(4.19) 



(m(p'J ® u{pi)) (l (g) 7^ - 7^ ® l) {u{p[) (g) u{pi)) = tiip' + m){p + m)'j^ 



p' x p 
PoPo 

Hence, at three level we find 

{k, k[\S\p,p[) = ^ {p' x p) S.{p^,p[, k,, k[) . (4.20) 

Including the non-scattering part ()4.14j) . we can read off the S-matrix as 

Sip,p') = l + ^(p'xp) + 0{g') , (4.21) 

or in terms of rapidities ()4.12|) as 

S{e, e') = l + 2ig smh{e - &) + 0(/) . (4.22) 

It is in fact possible to compute the S-matrix to all orders in g. For this computation 
to be possible, it is essential to quantize the theory in the pseudo- vacuum ()4.8|1 . This 
is because only in the pseudo-vacuum the higher loop diagrams are given by the sum 
of bubble diagrams, cf. Fig. ^2 The full computation is nevertheless rather involved 
and therefore has been relegated to App. IB. 31 The higher order corrections are given in 
()B.21|) and lead together with the non-scattering part to the exact S-matrix of the AAF 
model: 

Recall the value of the parameter g and m in terms of the 't Hooft coupling A from ()4.2|) . 

On the basis of integrability, the knowledge of the two-body S-matrix is enough to 
immediately write down the Bethe equations: 

giLmsinhe, _ TT 1 " sinh(% - 9k 



J-J-l + zosinhfØ, -^.) ■ ^ ' ' 



^^,. + igsmh{ej-ek) 

Energy and momentum are given by 

E = m ^ cosh Oj , P = m ^ sinh Oj . (4.25) 
i j 
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It is interesting to compare our result to the semiclassical spectrum of the AdS string. 
Because of the fermionic nature of the su(l|l) subsector, there are no classical spinning- 
string Solutions, but one can study states with very low momentum p ~ l/L (BMN 
states), which correspond to the lowest string modes. l/L corrections to the energies of 
these states were calculated exphcitly in fTB| and can be encoded in a compact form in 
a set of Bethe equations with the S-matrix^ 



S{p,p') = exp 

where^ 



- {p{e{p) - e{p)p') 



(4.26) 



A 

The same Bethe equations were shown to arise from quantization of the free fermions 
in the hght-cone gauge [40j. We now compare the S-matrix ()4.26p to ()4.23|) . For easier 
comparison, we give up our vector notation [p] and use p for the momentum and e for 
the energy. Furthermore we have to take into account that we considered the AAF model 
with rescaled time coordinate, i.e. we now compensate for this by rescahng the energy 
by so the dispersion relation becomes precisely ()4.27j) and the S-matrix is 

SKP.P) = ■ , • (4-28) 

1 - 4 \P^\P^) - ^{P)P ) 

It is immediately obvious that the S-matrices coincide in the low momentum approxi- 
mation, but deviate from each other at large momenta. 

We should warn the reader that the Bethe equations ()4.24|) cannot be used to describe 
the quantum spectrum of the string in AdSs x S^. Quantum corrections are obviously 
different for the quantum field theory defined by ()4.3|1 and the full string sigma-model. 
It is still interesting to note that quantum corrections in ()4.24|1 (deviations from the 
classical phase shift ()4.2fij) ) are not of the form anticipated for the quantum string in 
[ini 122] • We would also like to stress that ()4.24|) are bare Bethe equations, they have 
Solutions with negative energies, and the true ground state is the solution with the Dirac 
sea filled. 

We now turn to the discussion of the physical ground state of the AAF model. Recall 
that for technical reasons we quantized the theory in the pseudo vacuum |0), cf. ()4.8|1 . 
The physical vacuum \Q) is obtained from the pseudo- vacuum by exciting all negative 
modes, symbolically 

\n) = l[bHp,)\0)- (4.29) 

Pl 

This state corresponds to a non-trivial solution of the bare Bethe equations. Recall that 
the ^'''-particles are excitations with complex rapidities 6j = m — aj. Hence, the state 

^To compare to P0| one has to take into account that the R-charge denoted there by J is what we 
call L here. If one identifies L + (number of Bethe roots)/2 with the length, as is done in [201, then the 
S-matrix efFectively acquires an extra factor of exp[i(p' — p)/2], cf. eq. (4.31) in |2U| . 

^To be more precise, a lattice dispersion relation with replaced by 4sin^(p/2) was postulated in 
pn| . The difference, however, is 0(1/L^) for the BMN states. 
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Figure 9: Bethe roots for physical vacuum. With respect to the pseudo- vacuum |0), which 
underhes the "bare" Bethe equations, the physical vacuum \^}) is a highly excited state. It 
corresponds to an m-hne completely filled with Bethe roots aj. 



is described by placing Bethe roots as densely as possible along the ivr-line. This 
is depicted in Fig. El One can construct the vacuum solution of the Bethe equations in 
the thermodynamic (L —>■ oo) hmit, find the spectrum of physical excitations and their 
scattering S-matrix by solving certain integral equations. The details of the derivation 
can be found in the original literature on the massive Thirring model j42t 1^ , in the 
review or in the monograph ^^l- For our model, the density of roots in the vacuum 
Pn(tt) satisfies the foUowing equation 



mcosha = 2np^{a) 



+00 



da pn{a) 



2g cosh(Q; — a) 
I + g'^ sinh^ 



a — a] 



The energy density of the physical vacuum is given by 

E 

T 



— = — m 



da pn{a) cosh a 



(4.30) 



(4.31) 



The physical excitations are obtained from the vacuum configuration of Fig.IHlby inserting 
further Bethe roots outside the ivr-line or by making holes in the Dirac sea. The energy 
of those physical configurations is then measured with respect to the ground state energy 
dOH). 

Let us end this section with discussing the solutions of ()4.30|) in various regimes of the 
coupling g = In the attractive regime, g < 0, we find that pQ is determined only up 
to a constant. This is because a constant shift of the density: pn{a) — > pn{a)+C does not 
affect the left hand side of the equation. Choosing this constant arbitrarily large one can 
make the energy 1)4.311) arbitrarily negative, which leads to a vacuum instability for any 
attractive coupling constant. This behavior is similar to what happens in the massive 
Thirring model at infinite coupling, which corresponds there to the phase transition 
point. In the AAF model, however, the instability is present for all values of g < 0. 

In the repulsive regime, g > 0, there is no such problem. We solve ()4.3()j) by Fourier 
transformation. However, we are required to introduce a rapidity cut-off |a| < A in 
order to regularize the Fourier transform of cosh a. Then the root density of 1^2) can be 
written as the Fourier integral 



m 
16^ 



da e 



cosh I a 



cosh cosh 



g(l+ia)A ia)A' 
— + — 

1 + ta 1 — ta 



(4.32) 
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Figure 10: Coupling constant in AAF model. It is convenient to parameterize the coupling 
constant g hy a complex variable ^ in the way shown here. This is a solution of sin /i = l/g. 



Here /i = arcsin(l/5f), where we use the branches of arcsin as depicted in Fig. EU Notice 
that fi is real for strong coupling g > 1, but complex for weak coupling < g < 1. 
The integral ()4.32|) can be computed by contour integration. In the limit A — >■ cxd, the 
dominant contribution to the integral stems from the poles at a = in the strong 

coupling regime and from the poles at a = and a = ±i ^^^^ in the weak coupling 

regime. For g > 1 we have 

, , m A 1 TT^ TTa 

Pc a = 7^ e ^-/^ - cot — cosh , 4.33 

ZTT /i 2(/i — tt) n — fl 

for < (7 < 1 we have ()4.33|1 plus it complex conjugate. In the strong coupling regime, 
we can introduce a renormalized mass 

rriR me' '^-f^ (4.34) 

in order to get a finite result when the cut-off is removed. In the weak coupling regime 
this is not possible as the two terms in pn(a) have a different dependence on A. Hence, 
the model is non-renormalizable for < = < 1. It is instructive to express ()4.34|) 
in terms of the original coupling g. In the limit (7 — > 0, we find that the anomalous 
dimension of the mass m behaves as 

7m = ~ 7^ ■ (4.35) 

TT — jj, mg 

The approach of the anomalous dimension to zero as (? — > is not analytic. We see that 
the breakdown of renormalizability is a non-perturbative effect. 

In the strong repulsive regime {g > 1), the AAF model seems to be a well defined 
quantum theory. Unfortunately, additional complications arise in this case, because one 
can further diminish the vacuum energy by inserting n-strings around the ivr-line. We 
thus expect that the true vacuum is a condensate of n-strings, where n can depend on 
g, like in the strongly repulsive Thirring model 
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5 Faddeev-Reshetikhin model 



In this section we study string theory on x R. We should probably explain what 
we mean by that because S*^ x i? is not a string background. There are no problems 
with classical strings, in fact x R can be regarded as a subspace of AdS^ x S^, 
but quantization leads to UV divergences and non-zero beta-function. Althougli the 
resulting model cannot be interpreted as string theory, it is an interesting example of 
two-dimensional integrable field theory first considered in pU] . 
The string action in the conformal gauge is^ 



S = I drda 

An 



(5.1) 



where X° is the time coordinate and 

Ja = g-^dag . (5.2) 

Here g is a. group element of SU{2) that parameterizes an embedding of the string world- 
sheet in S^. The current satisfies the foUowing equations of motion: 

daf = 0, 

dajb - dbja + [ja,jb] = . (5.3) 

In addition we should impose the Virasoro constraints: 

ti jI = -2 {d^X'Y , (5.4) 

where d± = 9^ ± are derivatives with respect to the light-cone coordinates = 
(r ± (j)/2. Accordingly, j± are the light-cone components of the su{2) current. 

The standard way to proceed would be to quantize the sigma-model and then impose 
the Virasoro constraints in the weak sense - as subsidiary conditions on the physical 
states. This procedure is inconsistent for the model at hand, for the reasons explained 
above. It is also not likely to work for the full string sigma-model in AdSs x because 
of the well-known problems with the conformal gauge for the Green-Schwarz superstring. 
Another (in fact, no less standard) approach to string quantization is to fix the gauge 
completely, solve the Virasoro constraints (eliminate transverse degrees of freedom) , and 
only then quantize. We thus choose the temporal gauge: 

X° = KT , (5.5) 

such that the Virasoro constraints become 

ti jI = -2k^ , (5.6) 

and are solved by 

j± = ii^S± ■ a , (5.7) 



*The string tension of the AdS string is \/A/27r. 
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where S+ and S- are three-dimensional vectors of unit norm: S"^ 



equations of motion 



d+St + Ke'^''SiS'l = 0, 
d^Sl + Ke'^''S^S'l = 0. 



1. The satisfy the 



(5J 



Now we face an immediate problem: What Hamiltonian should be quantized? The 
original Hamiltonian, canonical conjugate to r in (jS.lj) is set to zero by the Virasoro 
constraints. The way out is to use another Hamiltonian and another Poisson structure, 
proposed in [20]. They are more natural from the point of view of integrability, and 
are potentially relevant in the AdS/CFT correspondence ØEl]. The action in the path 
integral for the evolution operator associated with the Hamiltonian of jHO] is 



S 



d x 



K 



(5.9) 



where C± are Wess-Zumino terms as defined in fj3.3p . This will be our starting point. It 
is easy to check that the equations of motion ()5.8jl follow from variation of this action. 

To develop perturbation theory for the FR model, we will perform the same change 
of variables as in sec. El 



Si + tSl 



2Sl 



Sl 



(5.10) 



Upon this change of variables, the Wess-Zumino term becomes a canonically normalized 
first-order action for </>_: 



S 



- Kv/(i-iø+i^)(i-iø-P) + + (i</>+r + \<i^-?) 



2k\ 



This action can be east into a very concise form if we combine 
component commuting spinor: 



Then ()5.1H) becomes a Dirac-like action^ 



(5.11) 



and (j)- into a two- 
(5.12) 



S 



The covariant derivative contains a field-dependent chemical potential: 

Dq = do — im — ig øø , Di = di . 



(5.13) 



(5.14) 



^To make power-counting possible we introduce the mass m = k and the coupling constant g = k/2. 



21 



This action describes a single charged particle and its anti-particle with the dispersion 
relations 

e = \/ jP' \ m?' — ra (particle) , (5.15) 

and 

E = a/p^ + W? + m (anti-particle). (5.16) 

This first of these equations is the BMN formula jini- The mass gap for the particles 
is offset to zero by the chemical potential. The energy of anti-particles is shifted in the 
opposite direction, so that anti-particles decouple at low energies and momenta. The 
low-energy effective theory can be constructed by defining big and small components of 
the spinor: 

^ = ^0, X=^ø. (5.17) 
The Lagrangian in ()5.13|) takes the form 

£ = ^<^*9o<^ + 2X*'9oX + ^V5*5iX + ^X*'9iVP + 2m|xP + ^?(^*V + X*V- Ixl') • (5-18) 
Integrating out x we arrive at 

L = i v*d,^ - ^ \dM^ - ^ [{v*d,^f + (9i^V)'] + . . . , (5.19) 

which is the same as ()3.20|) . This is another way to see that the Landau-Lifshitz model 
is the low-energy effective theory for strings on x R. 

Returning to the sigma-model, we can now quantize ()5.13j) in the "wrong" vacuum, 
in which all anti-particle states are empty. Then the S-matrix can be computed by 
summing the bubble diagrams. If the pole prescription is non-relativistic, the field- 
independent part of the chemical potential in ()5.14j) can be eliminated by a shift of the 
integration variable fc°, and we can use the relativistic dispersion relation instead of the 
BMN formulas ()5.15|) . ()5.16|1 . Calculation of the loop integral for the pp' pp' scattering 
(the details can be found in appendix IT^ gives: 

{f' + m)(S>{f + rn) + {f + m)®{f' + m) 
= 8m^sinh(g-gO • ^'-'^^ 
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The S-matrix is^° 

{kk'\s\pp') = ip„p'^(2„)H4k,k\p,p') 



E 

n=0 



+ m) ® {p + m) + {p + m) 



+ {27rf 5^^\k + k' ~p-p')j (out| (7° ® 1 + 1 ® 7° - 27^ ® 7^) 



8m2 sinh(^ - 9') 

x zc/ (7° ® 1 + 1 ® 7° - 27" (g) 7„) 



I in), 



(5.21) 



where by |m) and (out| we denote bi-spinors 



|m) = u{p) ® u{p) + u{p') ® u{p) , 
(out I = u{p) <S) u{p') + u{p') ® u{p) . 

Observing that 

|Z* + m = u{j))u{j)) , 
and taking into account the symmetry of the vertex, we can replace 

1, 



m] 



I in) (out I 



[f + m) ® [f + m) + [f + m) ® (j 
Then the S-matrix becomes 

'ig (out| (7° (g) 1 + 1 (g) 7° - 27" (g) 7„) I in) 



s{e,e') 



n=l 1- 



16m2 sinh(^ - 6') 



(5.22) 
(5.23) 

(5.24) 
(5.25) 



Plugging the exphcit form of the wave functions u{p) and u{p) into the matrix element 
of the vertex, we finally get: 



cosh^ 



coth ■ 



coth ■ 



(5.26) 



It is easy to check that in the hmit of small momenta the S-matrix of the LL model 
(I3.3()j) is recovered. 

The S-matrix of the FR model is not Lorentz invariant. This is not surprising because 
the Virasoro constraints explicitly break the Lorentz symmetry of the originally Lorentz- 
invariant chiral field. Because of the lack of Lorentz invariance the S-matrix is not a 



^°We change our conventions for the wave functions of external states. Here 

e-e/2 



u 



describes pseudo-particles with both positive and negative energies. The spinors are normahzed as 
u{p)u{p) — 2m. 
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function of 6 — 6' , which makes rapidity parameterization not the most convenient one. 
There is another parameterization of energies and momenta which is much more useful 
in this particular case: 

cosh 9 ~^ 



m — 1 

— = sinh 9 = — — 
m — \ 



sinh ^ = -3 -. (5.27) 



The S-matrix takes an extremely simple form in these variables: 



^(a;,x') = ^^^^. (5.28) 

^ ' ' x-x' ^2ig ^ ' 



We can now write down the Bethe equations: 



( 2imLxj \ YjXj-Xk + 2ig 



V a; ■ - 1 ; IJ: Xj -Xk- 2ig 

The states with Bethe roots in the interval —l<Xj<l carry negative energy. In 
quantum theory all of the negative-energy levels should be filled. The repulsive {g < 0) 
and the attractive {g > 0) cases are very different in this respect. In the repulsive case, 
the roots cannot form strings and the ground-state density in the thermodynamic limit 
is determined by the following integral equation 

(1^ = - 21^1 j ^ _ yy ^ • (5.30) 

The singularities at x = ±1 require regularization and will lead to renormalization of 
the bare parameters. It would be interesting to solve this equation and to perform the 
renormalization explicitly. 

The situation is more complicated in the attractive regime of g > 0. In this case, 
both positive- and negative-energy roots can form strings with Ax = 2ig. A simple 
calculation shows that the energy of the n-string with the centre of mass between — 1 
and 1 decreases with n (becomes more and more negative) up to n ~ l/2g. We thus 
expect that the ground state is a condensate of n-strings with very large n ~ l/g. This 
is in a qualitative agreement with the results of Faddeev and Reshetikhin for the lattice- 
regularized model [30]. The vacuum of the lattice model is a condensate of 25'-strings 
with S —>■ oc. We will not attempt to solve the Bethe equations ()5.29|1 here, but we 
expect that the renormalized solution describes the quantized principal chiral field j31j . 
as it is the case for the lattice model of pUj . 

6 Conclusions 

We hope that our calculations illustrate several features of the Bethe ansatz that can 
be useful in solving string theory in AdSs x S^. The consistency of the approach we 
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used fully relies on quantum integrability which reduces the many-particle problem to 
superposition of two-body interactions. This allowed us to avoid complexities associated 
with highly non-linear interactions in the sigma-models. Let us take the LL model as 
an example. Its Lagrangian contains infinitely many non-hnear terms, but we needed 
only the quartic vertices to derive the Bethe equations and thus to reconstruct the full 
spectrum. In other words we managed to solve the model by analyzing just the small 
fluctuations around the "north pole" n = (0, 0, 1). To illustrate this point, imagine that 
the target space of the LL sigma-model (the sphere 5*^) is deformed near the south pole, 
for instance that the constraint r?^ = 1 is replaced by n"^ = /(jt-s). If /(I —s) = l + 0{6^), 
the quartic vertices of perturbation theory around n = (0, 0, 1) are still the same and 
we will get the same two-body S-matrix. The complete non-perturbative spectrum of 
the deformed model, however, is totally different. The full spectrum is certainly not 
determined by the two-body S-matrix, because the deformation completely destroys the 
integrability and multi-particle interactions do not factorize any more. 

To further illustrate the power of integrability, we note that the Bethe ansatz com- 
pletely determines the spectrum of a system in a finite volume, which is particularly 
important for applications to string theory, but we did not use periodic boundary condi- 
tions anjrwhere in deriving Bethe equations. In fact, the S-matrix is only defined in the 
infinite volume where the notion of asymptotic states makes sense. 

Although Bethe ansatz completely determines the spectrum, it may happen that 
physically interesting states with low energies are described by rather complicated So- 
lutions of the Bethe equations. The ground state for a relativistic system typically is 
a Dirac sea that contains a non-trivial distribution of infinitely many roots^^. Usually 
one can find the ground-state distribution in a closed form only in the thermodynamic 
limit. It is interesting to note here that anti-particle-like states with negative energies 
seem to also arise on the gauge-theory side of the AdS/CFT correspondence. The all- 
loop Bethe equations for the spectrum of local operators in A/" = 4 SYM jH] do have 
negative-energy solutions at the non-perturbative level [2^1 EI- An interpretation and 
a physical significance of these states remains unclear to us. In particular, in J7j the 
negative-energy states were projected out rather than filled. We believe that these states 
are direct counterparts of the anti-particles in the sigma-model. Understanding their 
role in the AdS/CFT correspondence is a very interesting and important problem. 

Note added. A complementary approach to the quantum Bethe ansatz in the string- 
like sigma-models, which is based on the known exact physical S-matrices ^Hl, was 
developed in a parallel publication [Tfj . 
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A Computational details for LL model 
A.l Loop integrals 

A generic higher loop diagram in the LL model consists of a chain of bubbles as depicted 
in Fig. m In this appendix we compute the "bubble propagators" Ir- More specifically, 
we compute the momentum space representation of two parallel propagators D[t,x), cf. 
()3.16p . with two inflowing on-shell momenta p and p' {p > p'). We need to consider the 
cases with r = 0, 1,2 pairs of derivatives acting onto the two propagators. In the case 
without derivatives we have 




P 
p' 



p + p' - q 

{A.l] 

dtdx {D{x,t)f ei{p'+p")t-iiP+P')x 



du dq i 



2 



(27r)2 - q2 Jf- ^gj ^(^p2 p/2 _ ^) _ _^ _ g)2 _^ ^gj 

First perform the energy integral over ut by contour integration. The integrand has one 
pole in the upper and one pole in the lower half plane. Closing the contour in either half 
plane leads to 

, ,s i f dq 1 

h[p,p) = -7; TT 7 i r-K-2 w • 

2 J 2n [q- ^p- ^p'y - \{p - PY - le 

Shift the integration variable by g — >■ q + \p+\p' and perform the momentum integration 
also by contour integration. One finds 

h{p.p') = ^, ^ • (A.3) 
2{p- p') 

Loops with derivatives, see (jA.8j) and (jA.9|l . are similar to the above. Every pair of 
derivatives introduces a factor of —k{p — k) into the numerator. The energy integral in 
unchanged. However, the momentum integral is now divergent. For the evaluation of 
these integrals we invoke dimensional regularization. This has the same effect as adding 
appropriate counterterms to the action. The integrals become finite and their values are 
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A. 2 Summing all Feynman diagrams 

We compute 



{kk'\S\pp') 



{kk'\- 
n\ 



-|) jdxdt[{ip*d^ipf + {d,ip*ipY] 



\PP') ■ (A.6) 



The vertices are to be connected as in Fig. |3J Any other diagram will be zero. There are 
2"+^n! different ways of connecting n four-vertices among themselves and to the external 
legs in the way of Fig. ^ Furthermore we need to take into account that we have two 
different vertices. This amounts to placing derivatives onto the legs of the vertex, either 
onto the pair of in-going or on the pair of out-going lines. This leads to 2" different 
diagrams. 

A generic diagram consists of a chain of loops. If the derivatives hit the external 
legs we get a factor of —pp'. Derivatives on internal lines lead to three different kinds of 
loops: 



h{p,p) 
h{p,p') 







p 






P' 






p 


x: 




P' 






p 






p' 




p 



p 



(A.7) 
(A.8) 
(A.9) 



The corresponding loop integrals are evaluated in the previous subsection. Counting how 
many of the 2" diagrams lead to a given expression 

i-pp')*'ltip,pVt'iP,PVtiP,P) (A.io) 

is not completely trivial in general. This counting, however, becomes unnecessary as we 
find 

Hp.p') = {-pp')Hp.p') , h{p,p') = {-pp'fh{p.p')- (A.ll) 

This means that every diagram, no matter how the derivatives are distributed, will 
contribute ^ 



-ppTi^M) 



{2{p-p')) 



n-l ■ 



(A.12) 



Multiplying all factors together 



{k k'\S\pp') 



2"+^! ■ 2" ■ 



2 {^gy 



pp 



' 2 



-pp 



P — P' 



h{p,p',k,k') . 



-S+{p,p',k, k') 



(A.13) 



For n = 1 we recover the tree-level result ()3.27|) computed in the main text. 
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A. 3 Two-particle Green's function 

Here we compute 



(0| T 0{t, x)O{0, 0) 10) = ^! e-^*+^^- G{n, P) , (A.14) 

where O is defined in fj3.39p . For that we need to compute the loop diagrams Jq, h and 
h off-shell: 

. f duj dk 1 /•+°° dk 1 



(27r)2 {uj -k^)[Q-uj-{P- ky] iTii k^ + {P-ky-Q' 

(A.15) 

The answer is the same flA.3|l - ()A.5j) . where p, p' are defined as the roots of the denomi- 
nator in the integrand: 

In particular, 

/o(fi, P) = , ^ . (A.17) 

Summing the bubble diagrams as before we find: 

which gives ()3.40p . 

B Computational details for AAF model 
B.l Spinors 

In the AAF model we work with the foUowing two-component spinors: 

= [ cosr^(pO ) ' = [ -sinr^(p,) ) ^^'^^ 

and their Dirac conjugates u = «"''7° and v = f^"7°. The angle r]{pi) is defined through 
the relation cot2ri{pi) = ^. We use the following explicit representation of the Dirac 
matrices 

These spinors obey the equations of motion 

{p - m)u{pi) = , {p + m)v{-pi) = , (B.3) 
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where po = ^/pf + fn^- They satisfy the completeness relations 

u{pi)u{pi) = ^{f + m) , v{pi)v{pi) = ^{f-m) 

Zpo Zpo 



u{pi)u{pi) = -v{pi)v{pi) = — , u'^{pi)v{pi) = v\pi)u{pi) = — . (B.6) 



Pl 



Po 



Po 



(B.4) 



and the orthogonahty relations 

u\pi)u{pi) = v\pi)v{pi) = 1 , u\pi)v{pi) = v\pi)u{pi) = (B.5) 

and 



B.2 Loop integrals 

We compute the "bubble propagators" Ir in the AAF model. This is essentially a repeti- 
tion of the computation in the LL model (App. IA.1|) where now the fermion propagator 
()4.9|) is used. Defining P = p + p' we find: 



h{p,p) 




p + p' - q 



P 



P' 



D{q)® D{p + p - q) 



P XP 

-^1 ®7., 



(B.7) 



h{p,p') 




p + p' — q 



-P 
-P' 




p + p' - q 



P 
p' 



p x p 
4P2 



(B.8) 
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p + p' - q 



P 
p' 



p x p 

4P2 



-q x {p- q)Y D{q) (g) D{p + p - q) 



+ f P2 (/> ® 1 - 1 ® , 



+ m^P^ 1 



In the computation we made use of 



<fq 



(27r)2 (g2 + 



(fq 



(27r)2 (g2 + /\)2 4 



(B.9) 



(B.IO) 

(B.ll) 
(B.12) 



(27r)2 (g2 + /\)2 

for A E R. These formulas can be obtained from taking derivatives of the more general 
integral with numerator e*^ ". This integral itself can be evaluated by contour integration. 



B.3 Summing all Feynman diagrams 



We compute 
{kk'\S\pp') 



{kk'\ 



w. 



TH? 



\PP') ■ 
(B.13) 

For the reasons discussed in the main text only bubble diagrams contribute. As opposed 
to the LL model, here the chain of bubbles is made of two distinguished strands which 
connect one in-going particle with one out-going particle each. Le. there are two cases 
(Fig. ITT|l: p connected to k and p' connected to k' or p connected to k' and p' connected 
to k. The strands describe the flow of the spinor indices. At the vertex they hit onto the 
unit matrix or 7^. The two different possibilities can be tåken into account by writing 
the factor 

P=l® 7^-73®! , (ip)^ = lp (B.14) 

at each vertex in Fig. ^2 If there was only the first vertex and if we neglect the derivatives 
for a moment, then all diagrams at order (7" are given by 

{u{k[)®u{k^)-u{ki)®u{k[))¥(h{p,p')¥) u{pi)®u{p[)S+{pi,p[,ki,k[) . (B.15) 
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p 

p' 

p 

p' 



Figure 11: Higher loop diagrams in AAF model. 



Now we take into account the derivatives. The corresponding situation in the LL 
model was particularly lucky as any distribution of the derivatives led to precisely the 
same expression. Here, this is not the case and we are confronted with an rather large 
combinatorial problem. As all cases have to be handled separately anyway, it turns out 
the be more convenient to even consider all orders in g at once. The entire higher loop 
series is written in four terms: 

• no derivatives on external legs 



:B.16a) 



{u{ki) ®u{k[) -u{k[) ®u{ki)) Viou{pi) ® u{p[) S+{pi,p[, ki, k[) 
= {u{pi) ® u{p[) - u{p[) (g) u{pi)) Vio u{pi) ® u{p\) S-{pi,p[, kl, k[) 

• derivatives only on the in-going legs 

[f x p) {u{ki) (g) u{k[) - u{k[) (g) u{ki)) Vu u{pi) ® u{p\) 6+{pi,p[, ki, k[) 

= {p x P) {u{pi) (g u{pi) - u{p[) (g u{pi)) Vu u{pi) (g u{p\) S-{pi,p[, kl, k'i) 

(B.16b) 

• derivatives only on the out-going legs 

[k' x k) {u{ki) ® u{k[) + u{k[) ® u{ki)) Voou{pi) ® u{p[) 6+{pi,p[, ki, k[) 

= [p' x P) {u{pi) (g u{p\) + u{p\) (g u{pi)) Voo u{pi) (g u{p\) S^{pi,p[, kl, k[) 

(B.16c) 

• derivatives on the in- and out-going legs 

{p' x P) {k' x k) [u{ki) (g u{k[) + u{k[) (g u{ki)) Vqi u{pi) (g u{pi) 6+{pi,Pi, ki, k[) 

= [p' x py{u{pi) (g u{p\) + u{p\) (g u{pi)) Voi u{pi) (g u{p\) 6^{pi,p[, kl, k[) 

(B.16d) 

It remains to find expressions for the V^'s which sum all possible diagrams. The 
possible placements of the derivatives along the chain of bubbles are in one-to-one corre- 
spondence with binary sequences. One sequence of length n corresponds to one Feynman 
diagram at order g"'. Let "0" stand for derivatives "on the left" of the vertex (i.e. towards 
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out-going particles) and "1" for derivatives "on the right" of the vertex (i.e. towards 
in-going particles). For a diagram without derivatives on the external legs this sequence 
has to begin with a "1" and end with a "0", hence the notation Viq, etc. Now, this 
binary sequence translates into a sequence of bubble propagators Iq (no derivatives), Ji 
(one pair of derivatives) and I2 (two pairs of derivatives) according to the rules 

...00... ^ ...(-/i)... 

...01...^...(+/o)... (B.17) 

...10... ^ ...(-/a)... 

...ll...^...(+/i)... 

where we have chosen to include the sign from the second vertex (Recall that the vertex 
with the derivatives on the left comes with a minus sign, ()4.11|) ). The sign for the very 
first digit in the binary sequence has to be tåken into account separately. 

The next step is to generate all those sequences. The main idea is to think of an 
arbitrary binary sequence as a succession of constant subsequences, e.g. the sequence 
1111000100011111000 is made of 6 constant subsequences. It is clear, that all sequences 
are obtained from joining an arbitrary number of constant subsequences of arbitrary 
length. The constant pieces lead to 



P/± := J2 ^"(iP/i)" , (B.18) 



n=Q 



where the minus sign is used for a sequence of zeros and the plus sign for a sequence of 
ones. In ()B.18j) we included for each vertex a matrix P and one power of the coupling 
constant G = Between these constant subsequences we have to insert a factor GP/q 
when the sequence changes from "0" to "1" and a factor —GFI2 when the sequence 
changes from "1" to "0". Taking all this together we find 

00 

V^o ■=+Yl G'^''^¥/+(-P/2)P/_ (-P/0P/+P/2P/-) V , (B.19a) 

fc=0 

00 

Vu := + ^G2'=+^P/+(-P/2P/-P/oP/+)V , (B.lQb) 

fc=0 

00 

Voo := - ^ G^^^+^P/^ (-P/oP/+P/2P/-)V , (B.19c) 

k=0 

00 

Vbi := - ^ G^''^^¥I^{+FIo)FI+{-Fl2FI_FIoFI+fF . (B.19d) 

k=0 

The sums in ()B.18|) and ()B.19|) can be evaluated explicitly as functions of p, p' and m. 
It is computationally convenient to replace everywhere 

4^4 = (ip)'4(ip)' . (B.20) 

This is possible because every is sandwiched between two P's which satisfy the property 
()B.14|) . We refrain from printing the explicit K's as the formulas are rather lengthy and 
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not very illuminating. However, after inserting them into ()B.16jl . evaluating the spinor 
products and multiplying everything by the Jacobian PoPo/p' x V (cf- fl4.17|) ) we find the 
fairly compact result 

°° 2^ v' x v 

{kk'\S\pp') = -^/^^ \ S4p^,p[,h,k[) . (B.21) 

C Computational details in FR model 

C.l Loop integrals 

We need to compute 

f (Pk f + f-t + m ^ + m 

with the non-relativistic pole prescription. Performing the integral over first we get 

Loo = r" — Pik;p,p') 

27Ti 2{p + p'y{k-pi-ie){k-p[+te) ' ^ 

The denominator P{k;p,p') is a cubic polynomial in k which is symmetric in p and p': 

P{k] p, p) = 2E [^7° - (P - k)-f^ + m] ® (m - k-f^) 
- [E^ + k^-{P~ kf] 7° ® (m - k-1^) 
+ [E^ + e-{P- k/] [^7° - (P - kh^ + m] ® 7° 
-2E(fc2 + m2)70 ® 7° , (C.3) 

where E = pq + p^, P = pi + p[. Its most important property is that 

P{pi;p,p) = {f' + m)®{f + m), 

P{p[;p,p') = (p + m) + m) . (C.4) 

The integral in ()U.2|) diverges and requires regularization. Any reasonable regularization 
(for instance dimensional regularization) closes the contour of integration symmetrically 
in the lower and upper half-planes. Evaluating the residues at k = pi, k = p[ we get 
eq. ()5.2()j) in the main text. 
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